Ma- 

uht 





(NAS A -CR- 149661) 73ERM0DYNAMIC PROPERTIES N77-18938 

OF A HIGH PRESSURE SUBCRITICAL tJF6/He GAS 
VOLUME (IRRADIATED BY AN EXTERNAL SOURCE) 

(Florida Oniv.) 43 p HC A03/MF A01 CSCL 20M Onclas 

G3/77 16278 



MAR 1977 
RECEIVED 

NASA STI FACILITY .•?/ 
INPUT BRANCH '■* 


ENGINEERING AND INDUSTRIAL EXPERIMENT STATION 

College of Engineering 
University of Florida 


> 


Gainesville 


THERMODYNAMIC PROPERTIES OF 

A HIGH PRESSURE SUBCRITICAL 

UF /He GAS VOLUME 
6 

(IRRADIATED BY AN EXTERNAL SOURCE) 


BY 

DAVID E* STERRITT 
TENNESSEE VALLEY AUTHORITY 

GEORGE T. LALOS 
NAVAL SURFACE WEAPONS CENTER 

RICHARD T. SCHNEIDER 
UNIVERSITY OF FLORIDA 


WORK SUPPORTED BY 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
NGL 10-005-089 


DECEMBER 1976 


DEPARTMENT OF NUCLEAR ENGINEERING SCIENCES 
UNIVERSITY OF FLORIDA 



PREFACE 


This work was supported by the National Aeronautics and Space Administra- 
tion, Grant NGL 10-005-089. 

The authors wish to thank Dr. Karlheinz Thom, NASA Headquarters and Dr. 

H. "Bert" Helmick of Los Alamos Scientific Laboratory for advice and sugges- 
tions . 


ii 



ABSTRACT 

A computer simulation study concerning a compressed fissioning UF^ 

gas is presented. The compression is to be achieved by a ballistic piston 

compressor. Data on UF_ obtained with this compressor were incorporated 

in the simulation study. As a neutron source to create the fission 

events in the compressed gas, a fast burst reactor was considered. The 

15 2 

conclusion is that it takes a neutron f]ux in excess of 10 n/sec cm to 

produce measurable increases in pressure and temperature, while a flux in 
19 2 

excess of 10 n/cm sec would probably damage the compressor. 


TABLE OF CONTENTS 


PAGE NO. 

I. INTRODUCTION 1 

II, CONDITIONS FOR SIMULATION 5 

A. Experimental Constraints 5 

B. Treatment of the Test Gas as a Real Gas 5 

III. THE NUMERICAL SIMULATION 8 

A. Physical Data " 8 

B. Square Wave 8 

C. Analytic Energy Profile 9 

D. Temperature Anomilie (Limit of Model) 9 

E. Limit on Viscosity 9 

IV. RESULTS OF NUMERICAL SIMULATIONS 10 

V. CONCLUSIONS AND RECOMMENDATIONS 26 

APPENDIX 28 


iv 


I. INTRODUCTION 


Advanced reactor concepts like the gaseous core reactor, the nuclear piston 
engine, and the nuclear pumped laser, employ fissioning UF gas. For this 

D 

reason, the neutron physics of a gaseous critical assembly and the optical and 
thermodynamic properties of a fissioning UF gas are of interest. 

D 

One way to study fissioning UF in a subcritical configuration is to use a 

o 

ballistic piston compressor in connection with a pulsed high flux neutron source. 

The U.S. Naval Surface Weapons Center pioneered research Ln this field 
in the early sixties [1] ♦ Successful numerical analysis of the phenomena 
observed in ballistic compressors was achieved by Takeo [2] in 1965 with the 
development of a computer program which calculates gas leakage past the piston 
and heat losses from the test gas to the w*lls of the compressor. The University 
of Florida has developed the Ballistic Compressor Computer Code (BCCC) to aid 
in the experimental determination of thermodynamic parameters of gases at high 
temperatures and pressures from measurements of volume, temperature and pressure [3]. 

The BCCC has been modified to accept a variety of energy (fission) deposi- 
tion profiles such as those produced by the God.iva reactor. The resultant 
effects on the equation of state of the gas mixture and compressor are predicted 
before, during and after such energy deposition. 

A ballistic compressor shown schematically in Figure 1 consists of four 
main parts: the reservoir, the piston release section, the tube, and the high 

pressure section. 

The reservoir was designed for a maximum operating pressure of 136 atm and 
was statically tested to 200 atm. The reservoir is sufficiently large that the 
driver gas pressure remains essentially constant during the entire piston 
stroke. The tube is 3.89 m long and has a 5-cm bore. The 5-cm bore high 
pressure section is 17.78 cm long, has a wall thickness of 2.57 cm. and contains 
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diametrically opposite windows for absorption and emission measurements. 

DRIVER GAS 


t 



Fig. 1. Schematic diagram of ballistic piston compressor. 


The high pressure section is designed for pressures up to 5000 atm. Since 
the high temperatures generated in the high pressure section cause a slight 
vaporization of the inner walls, resulting in emission from impurities, the 
inner bore is plated with chromium. Figure 2 shows the piston and cup seals 
in detail. To minimize gas leakage around the piston, two cup seals are 
employed. The gases before and behind the piston expand these seals against 
the bore thus minimizing gas leakage past the piston. The steel piston body 
contains two phosphor bronze bearings that provide the surface upon which the 
piston moves and a molybdenum piston head to prevent ablation by the hot 
gas. The radial clearance between the bearing surface and the tube bore is 
50 ym. 



Figure 2. Schematic Diagram of Piston 
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In operation the chamber behind the plunger of the piston release section 
is pressurized so that the plunger moves to the forward (left) position and 
seals the reservoir from the tube. The tuhe is then filled with the test gas 
or gas mixture, and its pressure and temperature are adjusted. Next, the reser- 
voir is filled with driver gas to the pressure necessary to produce the desired 
maximum pressure in the test gas. The compressor is fired by releasing the 
pressure behind the plunger in the piston release section. The reservoir pres- 
sure acting on the front of the plunger moves it to its rear (right) position 
removing the seal between the reservoir and the tube. Reservoir gas rushes 
through the ports in the piston release 'section and impinges on the back of the 
piston driving it swiftly down the tube. The seals on the rear of the piston 
prevent all but a very small amount of reservoir gas from leaking into the tube 
and mixing with the test gas during a shot. Similarly, during the peak pressure 
part of the compression cycle, some test gas leaks across the piston and mixes 
with the reservoir gas. After the first compression cycle the piston oscillates 
back and forth until friction brings it to rest with equal gas pressure on its 
front and rear sides. The desired measurements are made during the peak pressure 
of the first compression stroke by the instrumentation located in the high 
pressure section. 

A rough approximation of the variation of test gas density and temperature 
with test gas pressure can be obtained from the polytropic relations 



where y is the polytropic exponent, the ratio of specific heats. Equating the 
work done on the piston by the expanding reservoir gas to the work done by the 
piston on the test gas being compressed, assuming that the reservoir pressure, 
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P^, is constant, and substituting Eq. (1) gives the following equation: 


/Pjn.iiA 

! ’i( r _ i) , 

In /», \\/\ 


(vr • 


which shows that P max does not depend on the length of the tube or piston 

diameter, but is a function of only. P , P , and Y. For (P /P )»1, 

o r max o 


Eg. (2) simplifies to 




which illustrates that the peak pressure generated in the test gas has a power 
law increase with reservoir pressure; i.e., a relatively low reservoir 
pressure generates a high test gas pressure, and that low y gases produce the 
highest peak test gas pressure for a given reservoir pressure. Typically, 
a reservoir pressure of 30 atm will generate a test gas pressure of 1000 atm 
in a monatomic gas. The ballistic pirton compressor was used before to measure 
the ratio of specific heats [4] (y) for UF^. 

In this report — based on the previous experience — a numerical simulation 
of such a ballistic piston compressor, being subjected to an intense neutron 
pulse during the time of maximum compression of the UF -gas, is presented. As 
a neutron source, the use of a bare fast burst reactor, such as Godiva IV was 
assumed. The expectation is that fission will occur in the UF C , and a fissioning 

D 

UF gas will be formed. It is intended to study this UF gas experimentally. 

6 6 

For safety reasons, the numerical simulation was necessary before the actual 


experiment is undertaken. 
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II. CONDITIONS FOR THE SIMULATION 

A. Experimental Constraints 

The reactor burst (~100 Usee) is much shorter than the 0.1 seconds it 
takes the compressor to reach maximum pressure and temperature. The objective 
is to initiate the reactor burst at the ~l/2 millisecond of high temperature 
ard pressure. 

After the compressor is readied for operation, the reactor is placed in 
in its supercritical operating condition awaiting a random neutron chain or an 
initiating burst from a neutron gun. The compressor is then fired and as 
the piston enters the high pressure test section a signal is produced in 
a magnetic pickup. This signal is used after suitable delay to initiate the 
neutron gun which initiates the reactor burst. 

Timing becomes critical here, since the compressor must be fired as soon 
as possible after the Godiva reactor awaits the burst of the neutron gun, 
least a random neutron chain result in an early burst from the Godiva reactor. 

Special consideration must be given to two aspects of this experiment. 
First, the increased reactivity from the rapidly compressed UF must be isolated 

o 

from the Godiva reactor. This can be effectively done by placing a thermal 
absorber between the moderating material for the compressor and the Godiva 
reactor. 

B. Treatment of the Test Gas as a Real Gas 

Nonideal processes to be considered in the simulation include gas leakage 
past the piston and heat loss from the test gas to the walls of the compressor. 

The Van der Waals* equation of state is: 

(p+ a/V 2 ) (V-b) = RT 

where a and b are the Van der Waals 1 constants. 
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the leakage of test gas past the piston is approximated by: 

dn/dz = — r (r *~r ) [-jp -Ip (r - r ) 2 + V] 

Ma o o p 8n 3x o p 

dn/dz - leak/ e rate 

p - density of the leakage gas 

Ma - molecular weight of the leaking gas 

r - inner radius of the barrel 
o 

r - radius of the piston 
P 

ri - viscosity of the leaking gas 

V - piston velocity 

3P 

— - pressur ; gradient alon B the piston 

Heat loss through the walls of the compressor and end plugs (one of which 

is the piston) is approximated by the cylindrical heat loss formula: 

3T ,3 2 T . 1 3T. 

3T ' ° ( ^2 + 7 SF> 

Boundary effects between the assumed homogeneous, turbulent test gas and 
the compressor walls are approximated by: 

Q = /T wall - T gas T plug - T gas , 

* R R plug ' 


where Q - heat lost during time step 

At - time step in seconds 
T wall - temperature of compressor walls 
T plug - temperature of end plug (piston) 
R, R plug - resistance to heat flow 


,-l 


min 


Airxk 

2ff0xh 


x 

K 

h 


- length of bore occupied by test gas 

- conductivity of the test gas 

- convective heat transfer coefficient of the test gas 
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. m 0.8 (r+C ) 3 0.3 
. K3 ,nVD, v 

k ' r V <-TT-> 

3 - Calibration factor 

n - Moles of test gas 

Vo - volume of te. c t gas 

D - bore diameter 

r - ideal gas constant 

C - constant volume specific heat 


A simple model is used to describe the fast burst reactor: 


where: 


dn _ p-B 
dt £ n 

p - reactivity 

B - delayed neutrons 

l - neutron lifetime 

n - neutron 


This approximation is meaningful only for sufficiently large reactivity, 
that is, all neutron sources except prompt neutrons may be neglected [4], 
The "shock wave" shutdown mechanism is approximated by: 

p ’ p 0 " “ T 

Where a is the negative of the temperature coefficient of reactivity and 
T is the increase in temperature above its initial value. The maximum 

power of the Godiva fast burst reactor is given by: [5] 


n = fcv~/2aK 

. p-B 

where w « — jj— 

K - is the reciprocal heat capacity. 

The energy released per unit is then: 

n * 4^7 sech^(wt/2) 
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This is the equation used in the BCCC to model the Godiva reactor. 

The full width at half maximum is given by: 
t = 3.524/W. 

A nore detailed description of the algorithums used in the BCCC is prer 
sented in reference [3]. 


III. THE NUMERICAL SIMULATION 


A. Physical Data 

The objective o/ the anticipated experiment is to produce a hot dense 

gas mixture of uranium hexafluoride in the field of a neutron flux of sufti- 

eient intensity to induce measurable changes in thermodynamic and optical 

properties of the compressed UF gas. 

6 

The initial temperature and pressure of the uranium hexafluoride is 
achieved in the ballistic compressor. In the base case, the initial conditions 
of the compressor are taken to be: 


Reservoir Pressure 

Uranium Hexafluoride Partial Pressure 
Helium Partial Pressure 
Piston Gap Calibration Constant 
Convection Calibration Constant 
Heat Transport Calibration Constant 
Piston Length Calibration Constant 


= 240 psig. 

= 0.10 atm. 

= 0.90 atm. 

= 0.01 cm. 

=0.1 dimensionless 
= 0.2 dimensionless 
= 25.0 cm. 


In the BCCC computer codes calibration factors are used which relate 


the actual gas leakage past the piston clearance to the computed leakage assum- 
ing a nominal clearance. The calibration factors obtained in the nonfission- 
ing experiments are assumed to be valid also for the fissioning experiment and 
are so used in this simulation. 


B. S quare Wave 

The initial modeling of the neutron flux is with a step insertion of 


fission energy: 
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E (t) = E^ (ergs/sec) < t < 

Where [t^, t^] is ohe interval during which is inserted. 

C. Analytic Energy Profile 

Hetrick [6] developed a theoretical model of the reactor burst more suitable 
for this analysis. As discussed previously, the flux or fluence as a function 
of time is given by: 

<|> - £w^/2aK sech 7 ' (wt/2) 

the energy produced during time At is: 

2 2 1 

E = At noapw /2aK sech (wt/2) (ergs/ see) 

where: a - energy per fission (ergs/ interaction) 

2 

<|> - neutron flux (n/cm sec) 

At - time interval (sec) 

255 

n - number of U atoms 

9 

a - absorption cross (cm ) 

2 

Let b i = a a £w /2aK 
1 a 

and bj = wt/2 

The range of values the parameters and b^ assume are: 
b^ - [10^, 10 17 ] ergs/sec 

b 2 - [10 3 , 10 5 ] sec" 1 

D. Temperature Anomilie (Limit of Model) 

Earlier investigations with the square wave fission energy profile demon- 
strated the limits of the model. Specifically, the high temperatures induced 
by the fission energy exceeded those temperatures over which the equation used 
to model the viscosity was fitted. The result was temperature singularities, 
followed by the leakage of all test gas past the piston. 

E. L imit on Viscosity 

Further investigations with the square wave and all cases with the analytic 
energy profile were calculated with the maximum viscosity limited to the value 
experienced at 2000°K. 
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IV. RESULTS OF NUMERICAL SIMULATIONS 


One objective of the simulation was to determine the sensitivity of the 
major physical variables on the results of the experiment. 

The maximum pressure is sensitive to the piston gap (i.e., leaking 
clearance between piston and bore wall) . For this reason, the sensitivity 
of the maximum pressure and temperature to changes in the piston gap is 
examined in Figures 3, A-l and A-2 (Figure A-l and A-2 are in the appendix), 
and summarized in Table I. 


TABLE I 




Maximum 

Density @ P max 

Maximum Temperature 

f * 

Figure 

Gap (cm) 

Pressure (atm) 

(#/cm 3 ) 

(°K) 

1 ; 

j : 

A-2 

.0050 

1600 

6.9-10 20 

1643 

i s 


.0075 

867 

3.2-10 20 

1474 

$ 

1 

A-l 

.0100 

423 

1.3-10 20 

1283 



.0125 

213 

5.9-10 19 

1096 

; 

3 

.0150 

117 

3.0-10 19 

933 

1 

( 1 

Figures 

3a, A-la and 

B-la have to be 

read as follows: At 

time t = 0 the piston 

* * 


is at rest about 380 cm away from the point of maximum compression (the test 
section end) near the breach of the pressure tube (see Figure 1) . After the 
plunger valve is opened at the time t = 0, the piston accelerates and comes 
to a standstill about 110 milliseconds later. The position value at this time 
(the minimum in the position-curve) divided by the position value at time 
t » 0 ig the compression ratio. After t = 110 milliseconds the piston starts 
to accelerate again in the reverse direction (now towards higher position-values) 
and undergoes a recompression around t = 220, 
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T I fi E - li ILL I SECONDS 

Piston Position, Gas Temperature and Pressure as a 
Function of Time 


ORIGINAL PAGE IS 
OF POOR QUALITY 


12 


From the time t = 0 to t = 110 milliseconds — the time of maximum compression — 
the temperature starts at ambient (t = 0, T = 300°K) and increases to its maximum 
value (in case of Figure 3a T = 933°K) which coincides with the minimum position 
value. Later, when the piston recedes, the temperature drops accordingly and 
rises again at the second recompression. 

Figure 3b (also A-lb and A-2b) shows the concentration of He and UF^ in front 

D 

of the piston. If there were no piston gap the number of moles in front of the 

> 

piston should not change. However, due to the existing gap, some of the helium 

which pushes the piston leaks naturally ahead of it thus increasing the number of 

moles in front of the piston (test section side) . The UF,. which is only on the 

6 

test section side of the piston, of course, does not leak over to the helium side, 

since at this side of the piston the pressure is much higher. Therefore, the 

UF C concentration stays fairly constant up to 100 milliseconds. Shortly after 
b 

this time (see Figure 3a) maximum compression occurs and the UF rushes past 

6 

the piston gap into the driver section of the compressor. The UF,. concentration 

b 

in the test section drops by a factor of 6. The same happens to the helium 
which leaked into the test section during the compression stroke. 

For this reason, the size of the gap has a rather dramatic influence on 
pressure and temperature. This is indicated in Figure 4. From previous experi- 
ences with UF -He mixtures [4] one can say that piston gaps between 50 and 75 
6 

micrometers are physically realizable. 

If the compressor is to be operated with fissionable UF gas, the most 

b 

significant variable next to the piston gap is the neutron flux or fluence to 
which the compressed gas is exposed. The fissioning events in the gas, due 
to this external neutron flux, give rise to energy deposited into the gass 
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ir Figure 3b 
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which gives : 


E ( t) = <|>0aNoN(t) At/ccsh (b 2 [t-to]) r 
= b^N ( t) At/cosh^ (b^ [t-to] ) , 


= <f>Noaa, 


where 


b^ = BCCC input - ergs/sec 

2 

4> = maximum flux n/cm sec 

-2 2 

0 = fission cross section, 577 x 10 cm 


a = 180 MeV/fission 

* 

235 

N ( t) = moles U 

At = duration of BCCC time step-sec 


b 2 = neutron pulse width parameter-sec 
No = Avagadro 1 s Number . 


This energy input relation was used in the compressor simulation. 

Table II shows typical results of this calculation. The piston gap used 
for these calculations was 100 micrometers. 

TABLE II 




i 
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Figure 5 shows the effect of deposition of fission energy into the gas. The 

curve labelled "Energy” shows the fission energy deposited into the gas as a 

15 2 

function of neutron flux. It takes a neutron flux in excess of 10 n/sec-cm 

to raise pressure and temperature measurably above the values which are prevailing 

20 

due to the mechanical compression of the gas. Of course, a flux of 10 , if 
it could be made available, would stress the test section to its limits. The 
nuclear energy input would, in this case, be more than 2 orders of magnitude 
higher than mechanical energy input. 

15 2 

Details of the compression for a flux of 10 n/cm sec can be seen in 
Figures 6a and 6b. The absolute values of pressure and temperature are slightly 
higher than in the case without neutron flux. 

Other than this, the details of the compression are identical to the non- 
fissioning case. Only if a substantial amount of nuclear energy is inserted 
(100 times the mechanical energy),* deviation in the compression detail becomes 
obvious. € See e.g.. Figures A-3,A-4 & A-5. The nuclear energy is inserted 
almost instantly (within 100 ysec). Temperature and pressure rises are steeper 
than in the mechanical case. The retrograde motion o* the piston starts earlier 


and the backswing is much more violent. 



Energy (joules) 


Figure 5 
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Concentration of UF^ and He as a Function of Time for a 
Flux of 10^ n/sec-cm^ 
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Another important aspect of the successful operation of the ballistic piston 
compressor is the precision of timing needed to obtain optimum temperature and 
pressure rise. The maximum compression lasts only about 0.5 millisecond. 

The neutron burst, which has a duration of about 100 p sec, has to be placed 
within the limits of the maximum compression time. 

Figures 7 and 8 show the consequences of the neutron burst if placed 
at the wrong time. In case of Figure 7, the neutron burst is 10 milliseconds 
early while, in Figure 8, it is 15 milliseconds late. In case the fission 
energy is inserted too early, the maximum pressure reached is, of course, lower 
than when placed correctly while an insertion too late creates a separate 
temperature maximum. Accordingly, in Figure 7 one sees a steep rise in tempera- 
ture when the nuclear energy insertion takes place; in Figure 8, one can see 
the temperature rise caused by the nuclear energy separately. A careful analysis 
of the situation shows that it is advantageous to insert the nuclear energy a 
short time — 200 microseconds — before the maximum pressure is about to be reached. 
Figure 9 shows the maximum pressure obtained, as a function injection time. 
Injection time is the time when the nuclear energy is inserted in respect to 
the time when the maximum compression — due to the kinetic energy of the piston — 
occurs. An injection time of "0" means that the nuclear energy is inserted at 
the same time when the maximum compression occurs, while an injection time of 
- 1 millisecond means that the nuclear energy was inserted 1 millisecond 

before maximum compression occurs. Figure 9 shows that the optimum injection 
time is - 200 microseconds. 

For safe operation of the device information concerning the highest neutron 
flux to which the ballistic compressor can be exposed is vital. A numerical 
simulation to find this highest possible flux was undertaken. The criterion 
is that if after the maximum compression, the piston returns to its original 
position, a neutron flux which caused this is no longer acceptable. In Figure 
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10 a case is shown where this flux is exceeded. The piston returns fast to 

its original position and overshoots it. The flux computed for this situation 
19 

was 4.8 x 10 n/cm sec. The result is satisfactory as far as safety consider- 
ations are concerned. It is impossible that such a high neutron flux would 
be made available inadvertently , and even then the damage to the reactor would 
be much more severe than the damage done to the ballistic piston compressor. 
Details on the gas mixture for this case can be found in Figure A-7 in the 
appendix. 

The maximum pressure and density achieved during compression may be 

improved not only by increasing the flux but also by increasing the partial 
235 

pressure of U F . At room temperature the partial pressure of UF is about 
6 o 

a tenth of an atmosphere. However , heating the compressor (heating tape) will 

raise the UF partial pressure above a tenth of an atmosphere. 

6 

235 

Table III itemizes the results of several variations in the U F^/He 
ratio. 




Table 

III 


Pressure 

(atm) 

Maximum Pressure 
(atm) 

Density 
@ P max 
#/cm 2 

Maximum Temperature 
°K 

.05 

368 


5 . 5 - 10 19 

1610 

.06 

384 


7 . 0 - 10 19 

1548 

.07 

399 


8.6-10 19 

1491 

.08 

415 


1.0 io 20 

1438 

.09 

429 


1.2 10 2 ° 

1389 

.10 

442 


1.3 10 2 ° 

1343 

.11 

454 


1.4 ] 0 2 ° 

1299 

.15 

464 


„* 

1115 


♦The test gas volume went to zero. 
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V. CONCLUSIONS AND RECOMMENDATIONS 

To maximize temperature or pressure, the insertion of the fission energy 

should precede the nonfissioning pressure maximum by about 200 microseconds. 

From the standpoint of undesirable side effects, error in the timing should 

only result in loss of data and no damage to equipment. 

16 2 

The maximum neutron flux should not exceed 10 n/cm sec and, if available, 
should be slowly adjusted upwards as experience provides data on the expected 
maximum compressor pressure. Care must be taken in increasing the maximum 
neutron flux, as error could result in destruction of the compressor. 

Care in securing the piston's cup seals will help in attaining the highest 
compressor pressure practicable as well as helping to localize fission products 
ahead of the piston. 

16 2 

However, operating at a lower flux than 10 n/cm sec will produce measur- 
able pressure and temperature increases. 
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APPENDIX 


For explanation of figures, see main text. 
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